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Abstract 

Background: Microscopy image segmentation lays the foundation for shape analysis, motion tracking, and 
classification of biological objects. Despite its importance, automated segmentation remains challenging for several 
widely used non-fluorescence, interference-based microscopy imaging modalities. For example in differential 
interference contrast microscopy which plays an important role in modern bacterial cell biology. Therefore, new 
revolutions in the field require the development of tools, technologies and work-flows to extract and exploit 
information from interference-based imaging data so as to achieve new fundamental biological insights and 
understanding. 

Results: We have developed and evaluated a high-throughput image analysis and processing approach to detect 
and characterize bacterial cells and chemotaxis proteins. Its performance was evaluated using differential interference 
contrast and fluorescence microscopy images of Rhodobacter sphaeroides. 

Conclusions: Results demonstrate that the proposed approach provides a fast and robust method for detection and 
analysis of spatial relationship between bacterial cells and their chemotaxis proteins. 



Background 

Modern bacterial cell biology has been revolutionised with 
the use of fluorescent markers coupled with microscopy 
allowing the visualisation of sub-cellular localisation in 
the bacterial cell Generating contrast in these images to 
determine the cell boundary is achieved using a number 
of optical methods including Phase Contrast and Differ- 
ential Interference Contrast (DIC) microscopy both of 
which depend on light changing its properties as it passes 
through the sample. 

It is increasingly important to interpret microscopy 
images in a quantitative manner thus being able to recon- 
struct the cell boundary from these images is of great 
importance, allowing the location of fluorescent markers 
to be determined within the bacterial cell In particular 
this information can allow correlation of the position of 
proteins within the bacterial cell throughout the cell cycle. 

Phase contrast images provide a clear light dark bound- 
ary around the entire cell making it easy to determine 
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by computational methods [1]. Software has been devel- 
oped to analyse phase contrast fluorescent microscopy 
images [2] however this software does not enable the anal- 
ysis of DIC images. DIC microscopy has the advantage of 
no phase ring in the objective and a larger depth of field 
making it advantageous for fluorescence microscopy in 
bacterial cell biology where the signals are often from a 
small number of fluorphores. It also enables better time 
courses to be taken with less photobleaching with less illu- 
mination needed per excitation. A lower excitation time 
per image results in a reduced level of photobleaching 
which is dependent on the intensity and duration of the 
excitation light. However in DIC image the cell outline is 
much more difficult to extract being made of both light 
and dark regions. In this paper we describe a new method 
for reconstructing the cell boundary from a DIC image. 
This then allows quantisation of the fluorescence image. 

DIC imaging 

The width of the emission spectrum of a common flu- 
orophore allows only for a limited number of spec- 
trally distinct fluorescent markers in the visible spectrum, 
which is also the regime where CCD-cameras are used in 
microscopy. For imaging of cells or tissues, it is necessary 
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to obtain an image from which the morphology of the 
whole cell can be extracted. This can be achieved by 
differential interference contrast microscopy [1]. How- 
ever, typically, DIC images appear with a bas-relief profile 
caused by the gradient of the optical path length and a 
phase shift between the two beams. For the human eye 
these DIC images are easy to interpret, see Figure 1. How- 
ever, automatic image analysis of DIC scans with hundreds 
of cells of different shapes and partially weakly identi- 
fiable contours is difficult. In the direction of the DIC 
shear vector, the intensity distribution of the cell is char- 
acterized by a transition of bright to dark, resulting in 
a well-defined contrast. But in the perpendicular direc- 
tion to the DIC shear there is no contrast against the 
background, and hence a lack of information about the 
complete cell boundary. 

DIC image segmentation 

Standard image processing methods for DIC image seg- 
mentation such as thresholding or edge detection produce 
insufficient results [3-5], with discontinuous regions or 
edges and can only be used under additional assumptions 
for the shape or by template constructions [6]. More- 
over, the use of deformable templates has been proposed 
[7,8]. These are modelled closed curves, which are fitted 
to object boundaries in iterative processes. More recently, 
[9] has demonstrated that the combined analysis of local 
image entropy and local illumination intensity could be 
suitable to identify individual cells with high sensitiv- 
ity and specificity at low computational cost. However, 
all of these techniques are either limited to cell types 
with relatively constant sizes and shapes, or require rel- 
atively long processing times, making them unsuitable 
for high- throughput applications. The alternative to these 
techniques is to perform DIC image reconstruction and 
then apply much simpler and more robust segmentation 
techniques. 



DIC reconstruction 

[10] has presented an excellent description, discussion, 
comparison and critical evaluation of most existing DIC 
image reconstruction methods. Most common iterative 
solutions are based on line integration and deconvolu- 
tion [11], variance filtering and directional integration 
using iterative energy minimization [12], and rotational 
diversity [13]. The latter technique involves taking several 
rotated DIC images and combining them using itera- 
tive deconvolution. Non-iterative methods include direct 
deconvolution [14] and the half-plane Hilbert transform 
[15] which is a qualitative Fourier-space approach to inte- 
grating the phase gradient. Moreover, a combination of 
the Hilbert transform and line integration method has 
been explored by [10]. 

Here we propose a high-throughput bioimage informat- 
ics approach to detect and characterize complex bacterial 
cells and chemotaxis proteins. The developed approach is 
based on a DIC shear orientation detection, followed by 
DIC image reconstruction and local segmentation of bac- 
terial cells and chemotaxis proteins. Finally, an analysis of 
spatial relation between bacterial cells and their chemo- 
taxis proteins is performed. R. sphaeroides was chosen 
because this is a small bacterial cell, generally 2 [ \im\ long 
and has chemosensory proteins known to localise in both 
the cytoplasm and the membrane. It has also previously 
been shown that there is cell cycle dependent position- 
ing of chemotaxis [16] proteins making this an ideal test 
for this software. The small size mean that automated 
analysis is important to reduce potential error in human 
measurement. 

Results 

Validation of the cell segmentation procedures 
DIC shear direction estimation 

We tested the performance of the proposed approach for 
DIC shear direction estimation on images of synthetic 
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Figure 2 DIC shear direction estimation procedure. DIC shear 
direction estimation from image of spherical bead of diameter 
1 [fimlAn estimated orientation = 43.52 [deg] (red line) was 
calculated with a rotated linear structuring element defined by / = 20 
[pixels] and n = 1 80. Scale bar corresponds to 1 [ fim]. 



spherical beads. An illustration of such an analysis is 
demonstrated in Figure 2. 

DIC image reconstruction 

We have applied the DIC reconstruction approach in a chal- 
lenging bioimaging application, based on the description 



in Methods Section, to extract and analyse bacterial 
cells and chemotaxis proteins in images of Rhodobacter 
sphaeroides. 

To quantify the performance of the DIC image recon- 
struction approach, seven complex images were selected. 
Manual segmentation of 290 bacterial cells was performed 
by an expert using the Pencil tool in Adobe Photoshop 
[17], as ground truth {Cgt)- Cgt cells and cells detected 
by the proposed approach (Ca) were compared using the 
Dice coefficient [18]. In our case, the Dice similarity coef- 
ficient is defined to find the overlapped regions between 
Cgt and Ca cell bodies in an image. The Dice coefficient 
is given by: 



D : 



2\c GT nc A \ 

\C GT \ + \C A \ 



(1) 



D lies between 0 to 1. If C G t = C A , then D = 1.0 
(perfect overlap), and if Cgt does not overlap with Ca, 
then D = 0 (no overlap). Comparison of overall accuracy 
for all analysed images is presented in Figure 3. The Figure 
shows, for every image, an average Dice coefficient value 
and error bar indicating its standard deviation. For all 
290 analysed cells, the Dice coefficient average value was 
D = 0.8635 ± 0.0513, while values of D = 0.8043 ± 0.1096 
were obtained when comparing two manual segmenta- 
tions performed at different times by the same expert of 
the same image set. 

In addition to the Rhodobacter sphaeroides, a set of 30 
polystyrene spherical beads (PolySciences Inc) was used 
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Figure 3 Validation - dice coefficient. Comparison of computational vs. manual method for 7 images. Bar plots show average Dice coefficients 
and error bars indicate one standard deviation. 
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to validate the performance of the DIC image recon- 
struction approach. Each synthetic bead had a diameter 
of 1 [fim] and an area of 0.7854 [fim 2 ]. For all recon- 
structed beads, the area average value was 0.8007 ± 
0.0312 [fim 2 ]. 

Furthermore, to compare the relationships between 
actual bacterial cell body and the same cell body deter- 
mined with a use of the DIC microscopy, a set of 20 
bacterial cells have been stained and then imaged using 
fluorescence and DIC microscopy, see Figure 4. The 
stained bacterial cells, observed on fluorescence images, 
were segmented using the standard thresholding method 
described in the Image segmentation Section. The same 
bacterial cells, observed on DIC images, were recon- 
structed and segmented using the automated DIC image 
reconstruction approach. For all analysed cells, the area 
average values were 1.4748 ± 0.2785 [ fim 2 ] and 1.4737 ± 
0.2396 [ fim 2 ], for bacterial cells on fluorescence and DIC 
images respectively. 



Spatial relations 

The proposed procedure for the analysis of the spatial 
relations between the cell and its chemotaxis proteins, 
as described in the Image analysis Section, has been 
also reviewed on the same set of test images containing 
290 cells. To quantify the performance of the proposed 
approach for analysis of the spatial relations between the 
cell and its chemotaxis proteins, estimated centroids and 
centrelines, for Cqj and Ca cells, were compared using 
the distance error measure proposed in [19]. In case of 
a centreline, the distance error is defined as the average 
distance between each point on the centreline of Cqt and 
the corresponding closest point on the centreline of C^. 
Distance error analysis of the cell centroids and centre- 
lines, for all analysed images, is presented in Figures 5(a) 
and 5(b). For all 290 analysed cells, the distance error 
average values were = 0.0551 =b 0.0192 [fim] and 
e d = 0.0818 zb 0.0402 [ fim], for centroids and centrelines 
respectively. 
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Figure 5 Validation: distance error. Comparison of computational vs. manual method for 7 images. Bar plots show average distance errors 6j and 
error bars indicate one standard deviation, for centroids (a) and centrelines (b) respectively. 
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Additionally, obtained results of the analysis of the spa- 
tial relations between the cell and its chemotaxis proteins, 
for a set of selected cells, is presented in Figure 6. 

Discussion 

An image processing and analysis concept for general 
differential interference contrast microscopy image seg- 
mentation has been developed (see DIC images Section), 
based on the DIC image reconstruction method using the 
2D Hilbert transform with a direction of reference in the 
Fourier domain. 

To evaluate the performance of the proposed approach, 
we applied it to images of the small bacterium Rhodobacter 
sphaeroides with fluorescently tagged chemosensory pro- 
teins, and the results are presented in Results Section. 



Visual inspection of the results confirm the robustness of 
the proposed approach for bacterial cells extraction and 
analysis of spatial relationship between bacterial cells and 
their chemotaxis proteins (see Figures 7 and 6). Further- 
more, we compared quantitatively the results obtained 
by the proposed approach with the ground truth results, 
delineated manually by an expert. As shown in Figure 3, 
the performed evaluation, based on the Dices similar- 
ity coefficient measure, demonstrates the accuracy and 
effectiveness of the proposed approach. Additionally, to 
address an issue of the DIC imaging producing a relatively 
large boundary of an imaged object, the proposed DIC 
image reconstruction approach was validated by com- 
paring our values with ground truth measurements in 
synthetic beads. 



W _ 




0 




(a) (b) (c) 











(d) 



(e) 



(0 




(g) 

Figure 6 Results of the proposed approach for the analysis of the spatial relations between the cell and its chemotaxis proteins (a-h). 

Legend: (green) - cell boundary, (yellow) - chemotaxis protein boundary, (yellow circle) - cell centroid position (r), (yellow star) - chemotaxis protein 
position (q), (magenta star) - point on the centreline p Cj which correspond to q, (green cross) - arbitrary end of the cell. Scale bars correspond to 
1 [fim]. 
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(a) (b) (c) (d) 




(e) (f) 

Figure 7 Bacterial cell and chemotaxis protein body segmentation. Bacterial cell and chemotaxis protein body segmentation workflow: (a) DIC 
input image, (b) DIC shear direction estimation: an obtained orientation = 48.20 [deg] (red line) calculated by use of a rotated linear structuring 
defined by / = 20 [pixels] and n = 1 80, (c) DIC image reconstruction, (d) an overlay of the segmented cell body in red, (e) fluorescent input image, 
(f) an overlay of the segmented chemotaxis protein body in red. Scale bars correspond to 1 [/xm]. 



This work lays the groundwork for using DIC to pro- 
duce spatiotemoral maps of proteins within the bacterial 
cell cycle. The use of DIC gives the potential for work- 
ing with smaller levels of photobleaching and generating 
higher resolution maps for bacterial cell biology. This new 
method allows the application of high throughput analy- 
sis of low copy number bacterial proteins throughout the 
cell cycle. The multiple parameters measured allows the 
determination of how any movement or positioning varies 
with the cell cycle and age of cell. 

Future work will be focused on high-throughput mea- 
surements of the Rhodobacter sphaeroides observed on 
temporal sequences of images. We will also investigate 
the applicability of the proposed concept to detect other 
bacterial cells observed in DIC images. 

Materials 

Rhodobacter sphaeroides 

R. sphaeroides strains were grown aerobically in suc- 
cinate medium [20] at 30°C with shaking at 225 rpm. 
Mid log phase cells were immobilized on a thin layer 
of 0.8% agarose in succinate medium on microscope 
slides [21]. 

One strain JPA1558 (TlpT-YFP) was used [21] enabling 
the visualisation of the cytoplasmic chemotaxis clusters 
using fluorescent microscopy. TlpT is used as a marker 



for the cytoplasmic cluster and CheA2 for the membrane 
cluster. 

Imaging 

DIC microscopy and fluorescence images were acquired 
with a Nikon TE200 microscope and YFP filter set 
(Chroma, Rockingham, VT) and recorded with a cooled 
charge-coupled device camera (Hamamatsu Photonics, 
Hamamatsu City, Japan). A Nikon oil immersion lOOx 
objective was used with an ND of 1.49. The final reso- 
lution achieved is 0.065573^0.065573 [ \im\ per pixel DIC 
illumination was achieved using Nikon Eclipse TE200 DIC 
attachment. This uses two Wollaston prisms and two 
polarisers to generate DIC. The 2nd polariser is mounted 
in the emission filter wheel and thus is not present when a 
fluorescence image is being acquired. Each slide had up to 
three images taken and images were analysed from three 
independent days, to ensure any analysis methods was 
independent of slide position or any differences in slide 
preparation. 

Methods 

Here we introduce an image processing approach for 
high-throughput detection and characterization of com- 
plex bacterial cells and chemotaxis proteins as observed 
on DIC and fluorescence microscope images (see 
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Figure 1). The workflow of the proposed approach is pre- 
sented in Figure 8 and its execution is demonstrated in 
Figure 7. 

For the DIC images, first, a robust DIC shear vector 
orientation detection procedure based on morphological 
gradient with linear structuring element and directional 
statistics are applied. Estimated orientation is then used by 
a DIC image reconstruction method where a 2D Hilbert 
transform with a direction of reference in the Fourier 
domain is employed. Global thresholding followed by a 
morphological filtering of the reconstructed image allows 
the detection of the bacterial cells. Additionally, a review 
process for ensuring quality of cells segmentation proce- 
dure is implemented. 

For the fluorescence images, a combination of morpho- 
logical filtering and local thresholding allows detection of 
the chemotaxis proteins. 

Finally, an analysis of the spatial relationship between 
segmented cells and their chemotaxis proteins is 
performed. 



DIC images: bacterial cells segmentation 
DIC shear direction 

Let us consider an image /(p), where p = [#,j] r is a 
column vector representation of the spatial location. For 
the detection of object boundaries, we use the morpho- 
logical gradient of the image by a structuring element S 
[22], defined as the difference between the morphologi- 
cal dilation 8 and erosion e of the image by structuring 
element S: 



gs(I) = 8 S (I) - € S (I) 



(2) 
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Figure 8 The workflow approach. The workflow of the proposed 
bioimage informatics approach for bacterial cells and their 
chemotaxis proteins identification in differential interference contrast 
and fluorescence microscopy images. 



where the dilation 8 of an image / by a structuring element 
S is defined as a locus of points z such that S hits / when 
its origin coincides with z: 



8 s (I) = {z:S z ni^0} 
and the erosion e is defined as follows: 

€ S (I) = {z:S z Cl} 



(3) 



(4) 



Here we assume that the DIC shear direction, observed 
in DIC images (Figure 7(a)), can be estimated by the 
highest output of its morphological gradient g calculated 
using rotated linear structuring elements S^ >a ) , defined by 
length / and direction angle a. Therefore, for each angle a n 
defined as: 



180 • (n - 1) 

N ' 



V«G [0,7V] 



(5) 



where N is a number of directions, the morphological gra- 
dient of the image / is calculated and the sum of all its 
pixels is determined: 



pel 



(6) 



As can be noticed, a set {w an } has circular distribution 
properties, therefore, in order to estimate its maximum 
(which corresponds to its mean orientation a^), direc- 
tional statistics is applied. The estimation of is done by 
fitting a von Mises distribution function to the {w an } data 
using a multidimensional unconstrained non-linear min- 
imization method [23]. In directional statistics, the von 
Mises distribution is a continuous probability distribution 
on the circle. This distribution is a circular analogue of the 
normal distribution and is defined as: 



fM({^a n }W^K) 



e KCOs(a—a fX ) 

2tzBq{k) 



(7) 



where the parameters and 1 //c are analogous to fi and 
cr 2 (the mean and variance) in the normal distribution. Bo 
is the modified Bessel function of order zero. 

The performance of this procedure applied to DIC 
image in Figure 7(a) is demonstrated in Figures 9 and 7(b). 

DIC reconstruction 

For a given DIC image / with DIC shear direction given 
by the orientation a^, the DIC image reconstruction is 
computed by applying a directional Hilbert transform, 
as described in Multi-dimensional Hilbert transform 
Section. The DIC image is reconstructed as follows: 



I R = mint) 



(8) 



Obara etal. BMC Bioinformatics 201 3, 14:1 34 
http://www.biomedcentral.eom/1 471-21 05/1 4/1 34 



Page 8 of 1 3 



o gradient data 
von Mises distribution 




80 100 120 

Angle [deg] 

Figure 9 DIC orientation estimation, von Mises distribution fit to a gradient values computed from DIC image (Figure 7) by use of a rotated linear 
structuring defined by / = 20 [pixels] and n = 1 80. Estimated mean orientation is = 48.20 [deg]. 



where 7^. is a Hilbert transform of the image I and its 
direction of reference in the Fourier domain is given by 
the unit vector defined as: 



e = [cosia^sinia^)] 



(9) 



In order to reduce the side effects of the Hilbert trans- 
form, only positive part of the image Ir is considered: 



Ho* 



if/* > 0 
if/* < 0 



(10) 



The graphic representation of the above described pro- 
cedure applied to ID DIC-like signal, is presented in 
Figure 10. In this graph, the DIC-like signal (solid line) of 
a simple ID object is modelled. In such a case, the DIC 
signal has the same levels inside and outside the object. 
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Figure 10 DIC reconstruction approach. The proposed approach for DIC reconstruction applied to 1 D DIC-like signal. 
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Application of the Hilbert transform (dash line) to such a 
DIC signal, results in a signal which has a level increased 
inside and decreased outside of the object. A side effect of 
the Hilbert transform is the appearance of spurious signal 
minima on the external side of the object, which we elim- 
inate by keeping only the positive values (dotted line) of 
the Hilbert transform. 

Figure 7(c) shows the performance of the proposed DIC 
image reconstruction procedure. 

Image segmentation 

Application of the DIC image reconstruction allows the 
reconstructed image to be automatically analysed using 
global or local thresholding methods [24], which would 
not be applicable on the original DIC image. In our case, 
the image Ir is segmented by Otsus global thresholding 
method [25], see Figure 7(d). 

Image postprocessing 

Small elements in the segmented image It are removed 
by performing a morphological opening by reconstruction 
with a disk structuring element «S( r ) of radius r: 

I G = PS (r) (YS (r) (lT),lT) (11) 

where ys is an opening of an image / by a structuring ele- 
ment S defined as an erosion of / followed by a dilation 
with S: 

ys(I) = 8 S (I) 6j(7) (12) 

Ps(I,J) is defined as a morphological reconstruction of 
mask image I from marker image /, / c /, which is 
obtained by iterating geodesic dilation of / inside I until 
stability [22]. 

Review 

Touching cells can be separated using watersheds [26], 
gradient flow [27] or active surfaces [28] method. In our 
case, a Euclidean distance map, calculated from touching 
cell boundaries, is used to calculate the watershed lines 
[29] which separates touching objects. 

Fluorescence images: chemotaxis proteins segmentation 
Image preprocessing 

In order to enhance low-contrast chemotaxis proteins 
within a fluorescence image 7, the top-hat filtering r 
with a disk structuring element S( r ) of radius r is 
applied: 

lH = r Sir) (I) (13) 
Image segmentation 

In order to segment the filtered image Ih> a local thresh- 
olding method based on the mean of the local intensity 



distribution is applied. In our case, the local neighbour- 
hood is represented by a circular window defined by 
a radius of 7 pixels. Resulting segmented image Ik is 
demonstrated in Figure 7(f). 

Image postprocessing 

Small elements in the segmented image Ik are removed 
by performing a morphological opening by recon- 
struction with a disk structuring element «S( r ) of 
radius r: 

Io = PS (r) (YS (r) (lK),lK) (14) 

Image analysis 

In order to analyse spatial relationships between bacterial 
cell and its chemotaxis proteins, the cells centreline has to 
be estimated. 

Cell centreline detection 

The most common approaches for centreline detection 
are based on Euclidean distance transform [30,31], fire 
propagation [32], Voronoi diagrams [33], clustering [34] 
and hybrid methods [35,36]. 

In our case, an approach which combines the Euclidean 
distance transform and the shortest path algorithm [37] 
is employed, and its workflow is shown in Figure 11 and 
Figure 12. In this approach, for every segmented cell C in 
the segmented image Iq, the cell boundary c is extracted. 
Then, the Euclidean distance transform of the cell body is 
calculated to define a cost map M. Afterwards, for every 
point on the cell contour cu the set of points k located 
near the half of the contour, in a range defined by s, 
is determined. Then, the shortest path p Ci ,c k . between C{ 
point and every point in the set k is calculated. Finally, the 
longest path Tc of all shortest paths is chosen. This path 
represents the centreline of the cell C. 

Where, EuclideanDistance(C) calculates an Euclidean 
distance map of every pixel of C to the nearest pixel 
outside of C (background pixel). ShortestPathiM^c^c^.) 



c i+f -s = c /ci 

'*> C *+f s+j = c k j 

/ c i+f+.s = C A- 2 , S+1 




Figure 1 1 Cell centreline detection procedure. A schematic 
diagram showing the cell centreline detection procedure. 
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Algorithm 1 Centerline extraction procedure. 
Require: C, s 
Ensure: Tc 

c <(— Boundary(C) 

M EuclideanDistance(C) 

M <- abs(M - max(M)) 

n size(c) 

for i = 1 — >• n do 

Ci <r~ c(i) 

k <- [i + f - s,i + f + s] 
for j = 1 2s + 1 do 

kj <- 

if fcj > n then 

end if 

c fcj <- c(kj) 

s Pci,c k . <— ShortestPath(M, Ci, c^) 
if d < length(sp Ci)Ck .) then 
d <— length(sp CtyCk ) 

end if 
end for 
end for 



Figure 12 Algorithm: centerline extraction procedure. Centerline extraction procedure. 



calculates minimum cost path, using Dijkstra shortest 
path algorithm [37], on a graph representation of the 
pixels of C having edges weighted by M values. 

Measurements 

A digital curve (arc) or path Fc (centreline in our 
case) is represented by a sequence of n distinct pixels 
p^P\>—iPn-\- The length element As; between two con- 
secutive pixels centres of the digital curve Tc is As; — 
\pi+i —pi\. The complete length of the digital curve which 
represents an open path is equal to [38]: 

n-2 

L (r c ) = L (po,Pn-i) = J] \Pi+l ~Pi\ ( 15 ) 
i=0 



Conclusion 

Presented work lays the groundwork for using DIC to pro- 
duce spatiotemporal maps of proteins within the bacterial 
cell cycle. The use of DIC gives the potential for work- 
ing with smaller levels of photobleaching and generating 
higher resolution maps for bacterial cell biology. This new 
method allows the application of high throughput analy- 
sis of low copy number bacterial proteins throughout the 
cell cycle. The multiple parameters measured allows the 
determination of how any movement or positioning varies 
with the cell cycle and age of cell. 

Future work will be focused on high-throughput mea- 
surements of the Rhodobacter sphaeroides observed on 
temporal sequences of images. We will also investigate 



The following measurements were implemented to 
determine the position of intracellular proteins relative to 
cell cycle (see Figure 13). 

• Cell centroid r, 

• cell centreline length I(r c )> 

• cell area and sum of cell intensity values, 

• middle point on the cell centreline p m , 

• clusters centroids c;, 

• point on the centreline p Ci which correspond to c;, 

• distance form p Ci to the furthest end 
max (L(pQ, Pcj ) t L(p cp p n -i))f 

• distance form p Ci to the arbitrary end, 

• distance form p c . to p m , L{p c . lPm y 

• cluster area and sum of cluster intensity values, 

• cluster minor and major axis lengths. 




Figure 13 Cell spatial measurements procedure. A schematic 
diagram showing the spatial measurements taken for bacterial cells 
and their chemotaxis proteins. 
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the applicability of the proposed concept to detect other 
bacterial cells observed in DIC images. 

Appendix 

Multi-dimensional Hilbert transform 

The Hilbert transform of one-dimensional function/ (#) is 
defined as follows: 

oo 

H{f(x)}=f n (x) = - [ —dr=f(x)*— (16) 

71 J T — X 71X 

— oo 

which, in the Fourier domain, is given by: 



F n (u) — F(u) • /sign(w) 



(17) 



where F and Fu are the Fourier transforms corresponding 
to / and fa respectively, u represents frequency and * is 
the convolution operator. 

However, in order to use the Hilbert transform with 
multi-dimensional functions, a direction of reference in 
the Fourier domain has to be introduced [39,40]. Hence, 
for a direction given by a unit vector e, the correspondence 



between a function F and its Hilbert transform F-u is 
defined as: 



Fn^) =i 7 (u) -isign g (u) 



(18) 



where u is a frequency coordinate and the multi-dimensional 
sign function is defined as: 



sign^(u) = 



+1 if u r e > 0 
0 if u r e = 0 
-1 if u r e < 0 



(19) 



In spatial domain, this correspondence is defined by 
convolution between / and the inverse Fourier transform 
of /sign- (u): 



^- 1 {/sign-(u)} = — ^f e (x) 



and we can write 



/^(x)=/(x)*^$f«(x) 

ttx 1 e e 
where x is a spatial coordinate. 



(20) 



(21) 
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Figure 14 DIC bacterial cells image analysis toolbox. Screen-shots of the DIC bacterial cells image analysis Toolbox 1 .0. 
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Software 

The DIC Bacterial Cells Image Analysis Toolbox (DICbc) 
has been implemented in MATLAB environment, see 
Figure 14. A detailed software description is presented 
below. Features: 

• Image Formats: 3 channels images - .jpg, jpeg, .png, 
.tiff. 

• Pre-Processing Algorithms: Hilbert transform, 
top-hat. 

• Image Thresholding: global - Otsu, local - mean, 
median, midgrey, Niblack, Bernsen, Sauvola. 

• Image Analysis: see Image analysis Section. 

• Post-Processing Algorithms: opening by 
reconstruction. 

• Supported Platforms: UNIX/Linux, MS-Windows 
2000/XP/Vista, Macintosh (OS X 10.1.4 and higher). 

• Online help documentation and a test data sets. 

Requirements: 

• MATLAB 6.5, R13 or higher. 

• Image Processing and Bioinformatics Toolboxes of 
MATLAB. 

• 512MB RAM. 

Installation instructions for MATLAB version: 

1. Create a directory in which you would like to place 
the software (e.g., C:\\DICbc or -/DICbc). 

2. Move the zipped toolbox to this directory. 

3. Unzip the software. This will create a directory called 
DICbc 1.0 which contains the Matlab code. 

Running instructions: 

1. After following the steps of Installation Instructions 
open a session of MATLAB. 

2. In the MATLAB current directory path change to 
the path where you install the toolbox. 

3. To run DICbc type in the MATLAB command 
prompt: run. 

The software is available upon request to: bogus- 
law, obar a@durham. ac . uk 
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